function [sigma, theta, l] = test_test_gp()

sigmas = [0.01, 0.1, 1, 10, 100, 1000];
thetas = [0.01, 0.1, 1, 10, 100, 1000];
ls = [0.01, 0.1, 1, 10, 100, 1000];

best_mse = 10000

for s = sigmas
    for t = thetas
        for l = ls
            mse = test_gp(s, t, l);
            if (best_mse > mse)
                best_mse = mse 
                stl = [s, t, l]
            end 
        end
    end
end

